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FORWARD  MULTIPLE  SCATTERING  CORRECTIONS 
AS  FUNCTION  OF  DETECTOR  FIELD  OF  VIEW 


I.  INTRODUCTION 

Low  visibility  atmospheric  conditions  occur  in  the  presence  of  heavy 

concentrations  of  atmospheric  aerosols  such  as  dust,  smoke  and  fog. 

Scattering  of  laser  radiation  by  such  aerosol  clouds,  being  basically  a 

multiple  scattering  process,  is  very  difficult  to  predict.  A  fruitful 

approach  to  the  problem  has  been  a  Monte  Carlo  technique,  because  it  can 

handle  strange  geometries  as  well  as  inhomogeneities 1  3.  Recently,  considerable 

attention  has  been  paid  to  analytical  solutions  to  the  equation  of  transfer 
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in  a  form  appropriate  for  the  laser  beam  propagation  problem  .  Although 
exact  solutions  have  not  been  obtained  to  date,  there  are  some  special  cases 
where  simple  and  useful  approximate  solutions  to  the  equation  of  transfer  are 
available.  For  tenuous  distribution  of  scatterers,  the  first-order  multiple 
scattering  theory  can  be  used,  and  for  dense  distribution  the  diffusion 
approximation  is  appropriate.  If  the  particle  size  is  large  compared  with 
incident  wavelength,  the  energy  scattered  by  the  particle  is  largely  confined 
within  a  small  angle  in  the  forward  direction  and,  therefore,  by  employing  the 
small-angle  approximation  it  is  possible  to  simplify  the  equation  of  transfer. 

In  Ref.  8,  a  systematic  study  of  contributions  of  increasing  order  of  scat¬ 
tering  for  both  realistic  and  model  aerosols  has  been  conducted. 
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As  indicatad  by  Ishimaru^' tha  first-order  multiple  scattering 

approximation  is  applicable  when  the  density  of  scatterers  is  so  low  that 

the  diffuse  (incoherent)  intensity  is  considerably  smaller  than  the  reduced 

(coherent)  intensity.  This  will  certainly  be  the  case  if  optical  distance 

traversed  by  the  bean  is  much  smaller  than  unity.  However,  the  sane  weak 

fluctuation  case  is  also  encountered  in  tha  situations  where  the  receiver 

has  a  narrow  receiving  angle.  In  this  case,  the  amount  of  scattered 

intensity  entering  into  the  receiver  is  small  compared  with  the  direct 

coherent  intensity,  and  therefore  the  received  field  is  predominantly 

coherent.  The  effect  of  the  detector's  finite  field  of  view  on  the 

received  power  has  been  a  subject  of  comprehensive  investigations  related 

12-17 

tc  forward  scattering  corrections  for  optical  extinction  measurements 

The  purpose  of  this  paper  is  to  study  the  effect  of  a  finite  field  of 

view  on  the  intensity  and  the  received  power  of  a  laser  beam  undergoing 

multiple  scattering.  Our  analytic  approach  is  based  on  the  theory  of 
18  19-22 

Dolin  and  Fante  ,  summarized  in  Section  II.  In  Sections  II  and  IV, 

we  apply  this  theory  to  Gaussian  beams.  Numerical  results  relevant  to 
the  beam  propagation  in  a  water  cloud  and  model  aerosol  particles  are 
presented  in  Section  V. 

II.  THE  FANTE-DOLIN  THEORY 

Our  considerations  will  be  based  on  the  equation  of  radiative  trans¬ 
fer  for  the  radiance  (specific  intensity)  distribution  function,  1(0, r,z). 
In  the  small  angle  approximation,  1(0, r,z)  satisfies 

I  •  ||  ♦  ♦  01  •  aa  P<$  -  $•}  I(|', r,z)  d  J*  (!) 
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where  0  end  are  volume  extinction  and  scattering  coefficients  (ra-^)  for 

the  aerosol  medium;  r  is  the  component  of  the  position  vector  transverse  to 

the  z  axis;  and  0  is  the  transverse  component  of  the  unit  propagation  vector 

(Fig.  1  in  Refs.  8  and  9). 

13  i9„22 

In  the  work  of  Dolin  and  Fante  ,  which  applies  to  sharply  peaked 

*+■ 

phase  functions,  1(0', r,z)  in  the  integrand  on  the  right  hand  side  of  Eq.  (1) 
is  expanded  in  Taylor  series  about  0*  ■  <5.  More  precisely,  l(0,r,z)  is 
split  into 

1(0, r,z)  -  lloJ  (0,r,z)  +  V*1  (0,r,z)  (2) 


where  superscripts  o  and  s  refer  to  unscattered  and  scattered  radiance, 
respectively.  Then,  one  expands 

(s)  -►  (a) 

IV  '  (0,r,z)  -  I  (0,r,z) 

+  Z  (0*  -  <t»v  -  I<S)  (l.r,«) 

k  k  3*k 

4n  (f  -  0).  (0*  -  I<S>  ♦  •  •  •  (3) 

2  k  l  k  4  3\3,,1e. 


where  k  and  l  refer  to  the  Cartesian  components  of  the  0  vector.  Recall 

that  0  »  0  2  +  0  ?. 

x  y 

Equation  (1)  now  yields 


3l 


(o) 


3  z 


+  01 


(o) 


(4) 
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ON- AXIS  IN 


FOV,  rad 


Figure  1.  Normalized  intensity  on  the  beam  axis  as  function  of  detector 
FOV  corresponding  to  a  ■  46.80  rad”^,  Rp  ■  1.0  cm,  and  t  -  1.0. 
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( s ) 

while  I  is  determined  from  the  nonhomogeneous  equation 


*  3l(”  0ui  ,,2..  .  3t<sl 

0  -yy  --<*>%!  *  -JJ- 


+  a(l  -  to)  i 


(s) 


GUI 


P(<J>  -  $•)  il°'  (<D',r,z)  d  f  (S) 


where  a)  »  Gg/a  is  the  single  scattering  albedo;  and  <$  >  is  defined  as 


<$  >  - 


-*-2  2 

P<4>)  <t>  d  $ 


(6) 


In  deriving  Eq.  (S)  we  have  neglected  the  terms  of  higher  order  than  second 
in  Eq.  (3).  The  solutions  to  Eqs.  (4)  and  (5)  can  be  obtained  by  the  method 
of  characteristics  in  the  Fourier  space  of  the  variables  $  and  r.  With 
the  Fourier  transforms  defined  generically  as 


«  -  •* 


d2cfr  d2r  I($,r,z)  *  n’r> 


(7) 


nt) 


2  ■*  iZ'Z 

d  <J>  P(d)  e  ^  v 


(8) 


we  obtain  the  following  solutions 


*(o)  ^  «♦  -*■  -*■  — oz 

i  (£.-n,z)  -  i(£  +  nz,n,z  «  o)  e 


(9) 


and 


:(s) 


(£.n,z) 


*  ■+  ■*  ■* 

HZ  +  nz,n,z  *  o) 


Pit  *  H(z-z') J  e"°z 


Z 

exp  |  -  |  <$2>  \t  ♦  n  (z  -  z")  I2  +  G(l-O))  ]  dz"  j  dz' 


(10) 


for  the  unscattered  and  scattered  contributions,  respectively. 


III.  CASE  OF  COLLIMATED  GAUSSIAN  BEAM 

We  assume  now  that  in  the  single  scattering  theory  the  scattering  phase 
function  P(0)  is  given  by  a  Gaussian  function,  i.e., 


-*•  t  2  -*2 

P(0)  M  ~  exp  (-a  0  ) 


(11) 


and  that  the  incident  collimated  beam,  directed  along  the  z  axis,  has  a 

Gaussian  spatial  form,  i.e., 

—1  2  { 2  \  -*»  2  -+2 

1(0, r, 2  -  0)  -  Fqtt  y  0  (0)  exp  (-Y  r  )  (12) 

where  5^  Is  the  Dirac  delta  function. 

In  the  situation  modeled  by  Eqs.  (11)  and  (12),  we  obtain,  after  performing 
the  inverse  Fourier  transforms  of  Eqs.  (9)  and  (10) ,  the  explicit  formulas 


for 

i(°) 

(0 ,  r,  z) 

( s )  -►  -*■ 

and  I  ( 0 , r , z ) . 

They  read 

I(o) 

(0, r,z) 

2 

_  -1  -a z  -y 

*  F  tt  e  e 

o 

r  5 (2>  (0) 

(13) 

and 

I(s) 

_  _  -az 

F  Oiiie  r 

2  -1 

(0, r, z) 

=»  - =—  dz ' 

(4  A( Z ' )  C(Z')  -B  (z')l 

(2ttT 

A(z')r2  -  3(z')0*r  C(z’)02 


4A(z')C(z*)  -  B  (z') 


where  the  functions  A,  B,  and  C  are  defined  as 


A(z') 


1  +•  au>z' 


B(Z') 


2z'  awz'2 


c (z ' ) 


1  z'2  Oujz'3/3 

„  2  *  2 
4y  4a 


(15) 
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In  polar  coordinates  (9,b),  the  vector  $  can  be  expressed  by 

<>  -  0  cos  b 
x 

$  "9  sin  b  (16) 

y 


For  a  detector  having  the  field  of  view  (FOV)  half-angle,  9Q,  the 

received  intensity  is  obtained  from  Eq.  (14)  after  integration  over  the 

angle  b  in  the  range  (0,2tt)  and  over  9  in  the  range  (0,9  ).  The  received 
( s ) 

intensity  F  (9Q,r,az)  corresponding  to  the  scattered  beam  thus  becomes 


F(s) (9D>r,az) 


°D 

I(3) (9,r,az) 9  d9 

0 


(17) 


where 


(  a) 

I  (9, r,az) 


-oz 

F  0)  e  rdz 

o _ 

2n 


(4A(|)C(f) 


B2(f>) 


Xo  [ 


9r  B (x/Q) 


4A(x/0)C (x/a)  -  B  lx/o) 


-] 


exp  (- 


A (x/a) r2  ♦  C (x/a)  92 
4A(x/o)C(x/a)  -  B2(x/a) 


wx 

e  dx 


(18) 


Here  I  is  the  modified  Bessel  function  of  zeroth  order,  and  x  is  a  dimensionless 
o 

integration  variable. 

The  contribution  coming  from  the  unscattered  part  of  the  beam  is 


F1®’  (S  -  F  -Tl  .'y2'2 

D  O 


(19) 


It  is  thus  seen  that  the  numerical  results  can  be  extracted  from  this 
theory  in  a  rather  simple  manner  by  performing  merely  a  double  integration. 
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IV.  CASE  OF  GENERAL  GAUSSIAN  BEAM 

More  generally,  one  can  account  both  for  the  spatial  and  angular  divergence 
of  t'ne  laser  bean  assuming  a  Gaussian  law  of  the  form 


—  22-2  2  2  2  -*2 

I(<J>,r,z  =*  0)  -  F  B  y  *  exp<-  8  -  Y  r  )  (20) 

o 

for  the  incident  beam.  Here,  the  parameter  B  describes  the  angular  divergence. 
3y  following  exactly  the  same  line  of  reasoning  as  in  Section  III.  we 

( S ) 

obtain  the  following  contributions  F  (9D,r,az)  corresponding  to  the  scat¬ 
tered  beam,  viz. , 


F  Is) (S  yr, az) 


(B’&D,yr, oz)  9'  d8‘ 


(21) 


where 


(si  -02 

I  (9,Yr,az)  =*  F  u  az  e  /2ir 
o 


fLr  2 

+  |  4A((x/a)  za}  c{ (x/a)  za}  -  b  { (x/a)  zo 

Jo  L 


r 


+ 1 


9 (yr)  -  B ( (x/a)  zo)/y 


4a{ (x/a)  za}  c{(x/a)  za}  -  B  {(x/a)  za} 


+  exp 


(7r)2  A((x/a)  za)/Y2  +  C((x/a)  zaiB2 
4A {(x/a)  za}  c{ (x/a)  za)  -  B2{(x/a)  za} 


-*■  expUo  xa  z)  •  dx' 

and  I  is  the  modified  Bessel  function  of  zeroth  order, 
o 


(22) 
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For  the  unscattered  bean,  we  obtain 


F*°*  (8D,Yr.or)  -  2Fq  S2  exp  (-Or)  6  2/n 


8'  d8*  i  {2(Y/<T>  8_8*  (yr)  oz) 

O  D 


+■  exp  {-  (Yr)2  -  £(y/o) 2  (or)2  ♦  B^J  9^  8’2} 


(23) 


The  integrals  in  Eqs.  (22)  and  (23)  have  limits  0  and  1.  This  choice  of 
limits  anticipates  the  use  of  Gaussian  quadrature  in  numerical  computation. 


V.  RESULTS  AND  DISCUSSIONS 

The  numerical  values  were  obtained  for  tha  on-axis  intensities  as  a 
function  of  8D,  and  tha  plots  shown  in  Figs.  1-3  depict  the  unscattered 
(reduced)  intensity  (squares) ,  scattered  or  diffuse  intensity  [circles) ,  and 
total  intensity  (triangles)  on  the  beam  axis  for  optical  depth  t  *  1,  4,  and 

IQ,  respectively.  A  divergent  beam  with  the  parameters  B  "  2ir/(YX)  end 

-1  24 

Y  ■  1.0  cm  is  assumed  to  be  propagating  with  the  Deirraend j ian  water  cloud 

model  Cl.  For  X  *  0.4S  ym,  the  extinction  coefficient  is  obtained  from 

the  Mis  theory  computations  when  a  modified  gamma  distribution  is  taken 

from  water  cloud  particle  site  distribution.  One  also  obtained  a  •  46.80  red 

for  the  Gaussian  fit  to  the  phase  function.  For  a  coaxial  detector,  with  a 

diameter  R^  -  1.0  cm,  the  same  as  the  laser  beam  diameter,  the  intensity  F 

becomes  a  function  of  8Q  and  s.  We  normalise  the  intensities  by  dividing 
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ON-AXIS  IN 


FOVt  rad 


Figure  3.  Normalized  intensity  on  the  beam  axis  as  function  of  detector 
FOV  corresponding  to  o  -  46.80  rad-1,  Rp  -  1.0  cm,  and  t  -  10.0. 
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out  the  factor  F  (8  -  00 ,  z  •  0).  For  the  saica  of  clarity  in  graphical 

presentation ,  we  also  scale  the  normalized  intensities  by  the  multiplicative 


factor  a:<p(r). 

In  order  to  investigate  the  significance  of  the  optical  thickness,  we 
calculated  the  power  received  by  a  coaxial  detector  of  radius  «  1  an. 
If  F(0Q,r,T)  denotes  generically  the  bean  intensity,  then  the  received 
power  is  defined  as 


9  <v 


2ff 


P 


F(8q,s,T)  r  dr 


C50) 


The  received  power  is  also  scaled  by  the  multiplicative  factor  exp(T), 
and  the  power  received  at  z  *  0  is  divided  out.  The  numerical  values  were 
obtained  for  received  power  for  a  coaxial  detector  of  diameter  ■  1  cm, 
and  the  situation  is  depicted  in  Figs.  4-6  corresponding  to  the  same  parameters 
as  in  Figs.  1-4,  respectively.  As  the  FOV  increases,  both  the  intensity 
and  the  received  power  saturate  rapidly,  independently  of  the  optical  depth. 
This  corresponds  to  the  situation  of  an  open  detector.  The  contribution  of  the 
scattered  power  becomes  dominant  in  the  saturation  region  for  optical  depth 
of  the  order  of  10. 

Gaussian  phase  functions  as  given  by  Eq.  (11)  were  best  fitted  in  Ref.  8 
to  the  exact  Mie  phase  functions  for  monodisperse  aerosols  with  radii  in  the 
range  from  2.01  to  40.2  pm.  In  Figs.  7-9,  we  show  the  intensity  vs. 
detector's  FOV.  Zn  addition,  we  show  the  received  power  vs.  FOV  for 
T  “  4.0  (Fig.  10).  These  sets  of  results  were  computed  for  the  model 
particles  characterized  by  the  parameters  ■  289.33  rad"1-,  and 


RECEIVED  POWER 


i  i  uvj 


Figure  4.  Received  power  (arbitrary  unira)  nn  v. 

Of  detector  FOV  corr.epondin^to  .  .  46.80  r.d'1,”  "“(.".'“L"0" 
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Figure  5.  Received  power  (arbitrary  units)  on  the  beam 
of  detector  FOV  corresponding  to  a  -  46.80  rad-1»  R_. 

t  -  4.0.  15 


axis  as  function 
■  1.0  cm,  and 
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DEIRMENDJIAN  WATER  CLOUD  MODEL  Cl 

A~  0.45  nm,  cl  =  46.8  rad"1 
<7=  1.633*  IO~3cm“*, 

Rl  =  1.0  cm,  Rq=  1.0  cm 

CASE  1  :  RD  =  Rl 

OPTICAL  DEPTH  T=  IQO 

-e-B-  UNXATTERED 
•e— e-  SCATTERED 
■A-A-  TOTAL 


0000  0.002  0004  0006  0.008  0010  0012 

FOV,  rad 

H  -.ure  Received  power  (arbitrary  units)  on  the  bean,  axis  as  function 
of  detector  FOV  corresponding  to  a  -  46. SO  rad-1,  r  -  l.o  Cn>,  and 


ON-AXIS  INTENSITY 
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FOV,  rad 

Figure  7.  Normalized  Intensity  on  the  beam  axis  as  function  of  detector 
FOV  corresponding  to  a  -  289.33  rad  ,  RD  -  1.0  cm,  and  t  -  1.0. 
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ON-AXIS  INTENSITY 


LARGE  SIZE  MONODISPERSIONS 

k  =  0.425  nm,  a  =289.33 rad'1 
cr=5iOxlO"3  cm'1, 

Rl  =  1.0  cm,  RQ=  1.0  cm 


CASE  i : 


Rp  =  Rl 


OPTICAL  DEPTH  T=4.0 

■e-e-  UNSCATTERED 
-e-e-  SCATTERED 
•Ar-Ar  TOTAL 


0000  0002  0004  0006  0.008  0010  0012 


FOV,  rad 

Figure  R.  Normalized  intensity  on  the  beam  axis  as  function  of  dete 
FOV  corresponding  to  a  *  ?89.33  rad-3,  R^  *  1.0  cm,  and  r  -  4,0. 


detector 


23 


V 


> 

H 

c n 

S  io‘ 


cn 

x 

< 

§  io2 


LARGE  SIZE  MONOOSPERSIONS 

A  =  0.425  urn,  a=  289.33  rad'1 
cr= 5.0x0"^  cm*1, 

Rl  =  1.0  cm,  Rq=  1.0  cm 


CASE  1 : 


=  R, 


OPTICAL  DEPTH  T=IQO 

-&-B-  UNSCATTERED 
■o— e-  SCATTERED 


0000  0.002  0004  0006  0.008  0010  0012 


FOV,  rad 


Figure  9.  Normalized 
FOV  corresponding 


malized  Intensity  on  the  beam  axis  as  function  of  detector 
ponding  to  a  *  289.33  rad-',  «  1.0  cm,  and  t  *  10.0. 


RECEIVED 


RQ  ■  0.2  cm  and  2.0  cm.  Results  were  obtained  for  the  same  parameters 


-2  -1 

a  -  5.0  x  10  cm  at  1*0. 45  pm,  for  a  coaxial  detector  of  diameter 
Rq  ■  1  cm.  Note  that  the  phase  function  and  the  Gaussian  fit  to  the  forward 
lobe  of  the  Mie  phase  function  are  plotted  in  Fig.  3  of  Ref.  8. 

In  order  to  study  the  effect  of  the  detector's  radius  R^,  we  let  Rp 
assume  values  smaller  and  larger  than  the  beam  diameter  of  the  1.0  cm, 
viz, 

as  in  the  case  of  R,,  -  1.0  cm.  However,  for  the  sake  of  clarity  we  present 
only  the  results  for  T  ■  4.0  for  the  two  cases,  viz,  9^  ■  0.2  cm 
(Figs.  11-12)  and  Rp  -  2.0  cm  (Figs.  13-14).  Figures  11-12  depict  the 
corresponding  values  of  intensity  and  received  power  as  functions  of  FOV 
for  RQ  ■  0.2  cm,  whereas,  Figs.  13-14  depict  the  corresponding  values  for 
Rp  *  2.0  cm.  The  three  cases  (R^  *  1.0,  0.2  and  2.0  cm;  T  ■  4.0)  for 
large  particles  with  phase  function  parameter  a  *  289.33  cm-1  show  clearly 
the  contrast  in  the  behavior  of  the  intensity  and  the  received  power. 

The  Fante-Dolin  approximation  employed  in  this  paper  enables  one  to 
estimate  the  corrections  to  the  Bouguer-Beer  law  for  a  receiver  having  a 
finite  field  of  view,  as  for  example  in  Ref.  25. 


VI .  RECOMMENDAT IONS 

1.  It  should  be  noted  that  the  small-angle  approximation  is  valid  for 
highly  forward-peaked  phase  functions.  It  is,  therefore,  recommended  that 
approaches  be  developed  to  deal  with  the  case  of  MS  in  laser  beams  traversing 
small  size  aerosol  particles,  with  broadly-peaked  phase  functions. 
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LARGE  SIZE  MONODISPERSIONS 

A  =  0.425  pm,  a= 28933  rad'1 
cr=5.0xl0‘3  cm"1, 

RL  =  1.0  cm,  Rq=0.2  cm 

CASE  2-  Rd  <  Rl 

OPTICAL  DEPTH  T=  40 


UNSCATTERED 

SCATTERED 


0000  0.002  0004  0006  0008  0010  0012 

FOV,  rad 

Figure  11.  Normalized  Intensity  on  the  beam  axis  as  function  of  detector 
FOV  for  the  same  parameters  as  In  Fig.  R  except  -  0.2  cm. 
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LARGE  SIZE  MONODISPERSIONS 

A = 0.425  jim,  a  =289.33  rad~ 1 
(7=50x10" 3  cm"*, 

Rl  =  1.0  cm,  Rq=0.2  cm 

CASE  2  :  Rd  <  RL 

OPTICAL  DEPTH  T=  40 


UNSCATTERED 

SCATTERED 


0000  0.002  0004  0006  0.008  0010  0012 

FOV,  rad 

Figure  12.  Received  power  (arbitrary  units)  on  the  beam  axis  as  a 
function  of  detector  FOV,  for  the  same  parameters  as  in  Fig.  11. 

The  received  power  is  multiplied  by  exp  (r). 
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OOOO  0002  0004  0006  0006  0010  00(2 


FOV,  rad 

Figure  13.  Normalised  Intensity  on  the  beam  axis  as  function  of  detector 
FOV  for  the  same  parameters  as  In  Fig.  8  except  “  2.0  cm. 
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FOV,  rad 

Figure  1A.  Received  power  (arbitrary  units)  on  the  beans  axis  as  a 
function  of  detector  FOV,  for  the  same  parameters  as  in  Fig.  13. 
The  received  power  Is  multiplied  by  exp(x). 
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2.  It  is  recommended  that  experiments  and  numerical  confutations  be 
performed  to  (a)  determine  the  effects  of  the  detector  field  of  view  on 

both  the  extinction  and  backseat taring  function  measurements  and  (b)  investigate 
the  effects  of  beam  diameter  on  optical  extinction  measurements. 

3.  It  is  recommended  that  the  above  investigation  be  repeated  for  the 
case  of  backseattering,  which  is  of  great  importance  for  single-ended 
electro-optical  systems. 
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ATTN:  DRDAR-SCP 

ATTN:  DRDAR-SCS 

ATTN:  DRDAR-TDC  (Dr.  D.  Gyorog) 

ATTN:  DRDAR-TSS 
ATTN:  DRCPM-CAWS-AM 
Dover,  NJ  07B01 

US  Army  Armament  Research  and 
Development  Command 
ATTN:  DRDAR-TSE-0A  (Robert  Thresher) 

National  Space  Technology  Laboratories 
NSTL  Station,  MS  39529 

Requirements  and  Analysis  Office 
Foreign  Intelligence  and  Threat 
Projection  Division 
ATTN:  DRDAR-RA I -C 
Aberdeen  Proving  Ground,  MD  21010 

Commander 

ARRADCOM 

ATTN:  DRDAR-QAC-E 

Aberdeen  Proving  Ground,  MD  21010 


37 


Dl rector 

USA  Ballistic  Research  Laboratory 
ARRADCOM 

ATTN:  DROAR-BLB  I 

ATTN:  DRDAR-TS8-S  1 

Aberdeen  Proving  Ground,  MD  21005 

US  ARMY  ARMAMENT  MATERIEL  READINESS 
COMMAND 

Commander 

US  Army  Armament  Materiel  Readiness  Command 


ATTN:  ORSAR-ASN  1 
ATTN:  DRSAR-IRI-A  1 
ATTN:  DRSAR-LEP-L  1 
ATTN:  DRSAR-SF  1 


Rock  Island,  IL  61299 
Commander 

US  Army  Dugway  Proving  Ground 

ATTN:  Technical  Library  (Docu  Sect)  1 

Dugway,  UT  84022 

US  ARMY  TRAINING  &  DOCTRINE  COMMAND 
Commandant 

US  Army  Infantry  School 

ATTN:  CTDO,  CSD,  NBC  Branch  1 

Fort  Bennlng,  GA  31905 

Commandant 

US  Army  Missile  &  Munitions  Center 
and  School 

ATTN:  ATSK-CM  1 

Redstone  Arsenal,  AL  35809 

Commander 

US  Army  Logistics  Center 

ATTN:  ATCL-MG  1 

Fort  Lee,  VA  23801 

Commandant 

US  Army  Chemical  School 

ATTN:  ATZN-CM-C  1 

ATTN:  ATZN-CM-AD  2 

ATTN:  ATZN-CN-CDM  (Dr.  J.  Scully)  1 

Fort  McClellan,  AL  36205 

Commander 

USAAVNC 

ATTN:  ATZQ-D-MS  1 

Fort  Rucker,  AL  36362 


Commander 

USA  Combined  Arms  Center  and 
Fort  Leavenworth 
ATTN:  ATZL-CAM- I M 

Fort  Leavenworth,  KS  66027 

Commander 

US  Army  Infantry  Center 
ATTN:  AT  SH-CD-MS-C 

ATTN:  ATSH-CD-MS-F 

ATTN:  ATZB-DPT-PO-NBC 
Fort  Bennlng,  GA  31905 

Commander 

USA  Training  and  Doctrine  Command 
ATTN:  ATCD-N 

ATTN:  ATCD-TEC  (Dr.  M.  Paste:) 

ATTN:  ATCD-Z 

Fort  Monroe,  VA  23651 

Commander 

US  Army  Armor  Center 
ATTN:  ATZK-CD-MS 

ATTN:  ATZK-PPT-PO-C 

Fort  Knox,  KY  40121 

Commander 

US  Army  TRAD0C  System  Analysis  Activity 
ATTN:  ATAA-SL 

ATTN;  ATAA-TDB  (L.  Dominguez) 

White  Sands  Missile  Range,  NM  880U2 

Commander 

USA  Field  Artillery  School 
ATTN:  ATSF-GD-RA 

Ft.  Sill,  OK  73503 

Dl rector 

USA  Concepts  Analysis  Agency 
ATTN:  M0CA-SMC  (Hal  Hock) 

8120  Woodmont  Avenue 
Bethesda,  MD  20014 

Los  Alamos  National  Laboratory 
ATTN:  T-D0T,  MS  B279  (S.  Gerstl) 

Los  Alamos,  NM  87545 
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US  ARMf  TEST  &  EVALUATION  COMMAND 
Commander 

U$  Army  Test  i  Evaluation  Command 
ATTN:  DRSTF-CM-F 

ATTN:  DRSTE-CT-T 

ATTN:  DRSTE-AD-M  (Warren  Baity) 

Aberdeen  Proving  Ground,  MO  21005 

r ommander 

i ;  S  A  fPG 

ATTN:  STEEP-MM-IS 

ATTN:  STEEP-MT-OS  (CPT  Decker) 

Ft.  Huachuca,  AZ  85613 

d  er 

Dugway  Proving  Ground 

ATTN:  STEDP-MT  (Dr.  L.  Solomon) 

Cugway,  'JT  84022 

DEPARTMENT  OF  THE  NAVY 

Commander 

Naval  Researcb  Laboratory 

ATTN:  Code  5709  (Mr,  W.  E.  Howell) 

ATTN:  Code  6532  < Mr ,  Curclo) 
attn:  Code  6532  (Mr.  Trusty) 
attn:  Code  6530-2  (Mr.  Gordon  Stamm) 

*  "  T  s :  rode  8320  (Dr.  Lothar  Ruhnke) 
a T T n ;  Code  8  326  (Dr.  James  Fitzgerald) 
attn:  Cede  43202  (Dr.  Hermann  Gerber) 

4555  Overlook  Avenue,  SW 
Washington,  DC  20375 

"'riel,  Bureau  of  Medicine  S  Surgery 
department  c t  the  Navy 

attn;  MED  3C33 
Washington,  DC  20372 

Comma  n  der 

Naval  *  I r  Systems  Command 
A1TN:  Code  AIR-3010  (Dr.  H.  Rosenwasser) 

ATTN:  Code  AIR-5363  (D.  C.  Caldwell) 
Washington,  DC.  20361 

Comma nder 

Naval  Sea  Systems  Command 

ATT, 4;  SI  A-62T 1  3  (LCDR  Richard  Gilbert) 

ATTN:  ST  A-62Y2 1  (A.  Kanterman) 

attn;  S(*-6?V2I  (LCDR  W.  Major) 

Washington,  DC  20362 


Project  Manager 

Theatre  Nuclear  Warfare  Project  Office 
ATTN:  TN-09C 

Navy  Department 
1  Washington,  DC  20360 

I 

1  Institute  for  Defense  Analysis 

400  Army-Navy  Drive 
Arlington,  VA  22202 

Commander 

1  Naval  Surface  Weapons  Center 

t  Dahlgren  Laboratory 

ATTN:  DX-21 
ATTN:  Mr.  R.  L.  Hudson 

ATTN:  F- 56  (Mr.  Douglas  Marker) 

Dahlgren,  VA  22448 

1 

Commander 

Naval  Intelligence  Suport  Center 
ATTN:  Code  434  (H.  P.  St.Aubln) 

4  30  1  Su 1 1 1  and  Road 
Sultland,  MD  20390 

1  Commander 

1  Naval  Explosive  Ordnance  Disposal 

1  Technology  Center 

1  ATTN:  AC- 3 

1  Indian  Head,  MD  20640 

1  Of f I cer- I n-Ch arge 

Marine  Corps  Detachment 
Naval  Explosive  Ordnance  Disposal 
Technology  Center 
Indian  Head,  MD  20640 

1  Commander 

Naval  Air  Development  Center 

ATTN:  Code  2012  (Dr.  Robert  Helmbold) 

Warminster,  PA  18974 

1  Commander 

1  Naval  Weapons  Center 

ATTN:  Code  382  (L.  A.  Mathews) 

ATTN:  Code  3882  (Dr.  C.  E.  Dlnerman) 
ATTN:  Code  3918  (Dr.  Alex  Shlanta) 
China  Lake,  CA  93555 

1 

1  Commanding  Officer 

1  Naval  Weapons  Support  Center 

Applied  Sciences  Department 
ATTN:  Code  50C,  Bldg  190 

ATTN:  Code  502  (Carl  Lohkamp) 

Crane,  IN  47522 
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US  MARINE  CORPS 

Commanding  General 
Marine  Corps  Development  and 
Education  Command 

ATTN:  Fire  Power  Division,  009!  i 

Quantlco,  V A  22134 

DEPARTMENT  OF  THE  AIR  FORCE 

Department  of  the  Air  Force 
Headquarters  Foreign  Technology  Division 


ATTN:  TQTR  1 

Wright -Patter son  AFB,  OH  45433 

HQ  AFLC/LOWMM  1 

Wr 1 gh t-Pat terson  AFB,  OH  45433 

AFAMRL/TS 

ATTN:  COL  Johnson  1 

Wr I ght-Patterson  AFB,  OH  45433 

HQ  AFSC/SDZ  1 

ATTN:  CPT  D.  Rledlger 

Andrews  AFB,  MD  20334 

USAF  TAWC/THl  1 

Egl In  AFB,  FL  32542 

USAF  SC 

ATTN:  AD/YQ  (Dr.  A.  Vaslloff)  1 

ATTN:  AD/YQO  (MAJ  Owens)  1 

Eg  I  I n  AFB,  FL  32542 

AO/SRO  1 

Egl  In  AFB,  FL  32542 

Dr.  Charles  Arpke  1 

OSV  Field  Office 

P.0.  Box  1925 

Egl  In  AFB,  FL  32542 


Commander 

Hanscom  Air  Force  Base 

ATTN;  AFGL-POA  (Dr.  Frederick  Volz)  1 

Bedford,  MA  01731 

Headq  uarters 
Tactical  Air  Command 

ATTN;  DRP  1 

Langley  AFB,  VA  23665 


AFOSR/NE 

ATTN:  MAJ  H.  Wlnsor  1 

Bolling  AFB,  DC  20332 

OUTSIDE  AGENCIES 

Battel le,  Columbus  Laboratories 
ATTN:  TACTEC  1 

505  King  Avenue 
Columbus,  OH  43201 

Toxicology  Information  Center,  JH  652 
National  Research  Council  1 

2101  Constitution  Ave. ,  NW 
Washington,  DC  20418 

Dr.  W.  Michael  Farmer,  Assoc  Prof,  Physics 
University  of  tennessee  Space  Institute  1 

Tullahoma,  TN  37388 

ADDITIONAL  ADDRESSEES 

Office  of  Missile  Electronic  Warfare 
ATTN:  DELEW-M-T-AC  (Ms  Arthur)  1 

White  Sands  Missile  Range,  NM  88002 

US  Army  Mobility  Equipment  Research  and 
Development  Center 

ATTN:  DROME-RT  (Mr.  0.  F.  Kezer)  1 

Fort  Bel  voir,  VA  22060 


Dl  rector 

US  Night  Vision  and  EO  Laboratories 
ATTN:  DRSE  L-NV-V I  (Dr.  R.  G.  Buser)  1 

ATTN:  DRSEL-NV-V  I  (Mr.  R.  Bergemann)  1 

ATTN:  DELNV-VI  (Luanne  Obert)  1 

ATTN:  DELNV-L  (D.  N.  Spector)  1 

Fort  Belvolr,  VA  23651 


Commandant 

Academy  of  Health  Sciences,  US  Army 
ATTN:  HSHA-SSI  2 

Fort  Sam  Houston,  TX  78234 
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